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Abstract 

Different components of heritability, including genetic variance (V G ), are influ- 
enced by environmental conditions. Here, we assessed phenotypic responses of 
life-history traits to two different developmental conditions, temperature and 
food limitation. The former represents an environment that defines seasonal 
polyphenism in our study organism, the tropical butterfly Bicyclus anynana, 
whereas the latter represents a more unpredictable environment. We quantified 
heritabilities using restricted maximum likelihood (REML) procedures within 
an "Information Theoretical" framework in a full-sib design. Whereas develop- 
ment time, pupal mass, and resting metabolic rate showed no genotype- 
by-environment interaction for genetic variation, for thorax ratio and fat 
percentage the heritability increased under the cool temperature, dry season 
environment. Additionally, for fat percentage heritability estimates increased 
under food limitation. Hence, the traits most intimately related to polyphenism 
in B. anynana show the most environmental-specific heritabilities as well as 
some indication of cross-environmental genetic correlations. This may reflect a 
footprint of natural selection and our future research is aimed to uncover the 
genes and processes involved in this through studying season and condition- 
dependent gene expression. 



doi: 10.1002/ece3.718 



Introduction 

Organisms inhabiting heterogeneous and/or seasonal envi- 
ronments often show phenotypic plasticity in which a sin- 
gle genotype yields different phenotypes in response to 
biotic and/or abiotic aspects of the environment (Pigliucci 
2001). The adaptive value of phenotypic plasticity is trait 
specific; for some traits expressing different phenotypes 
under heterogeneous environments maximizes fitness, 



whereas for others maintaining the trait value over a 
range of conditions (phenotypic canalization) is a more 
beneficial mechanism (Stearns and Kawecki 1994). In 
order for phenotypic plasticity to evolve, genetic variation 
is required in the form of genotype-environment interac- 
tion (Via and Lande 1985). Furthermore, understanding 
the evolutionary potential of organisms requires quantify- 
ing the amount of genetic variability expressed for the 
traits of interest that are under selection. Therefore, to 
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understand the evolution of phenotypic plasticity and to 
assess the potential for future evolutionary change, it is 
essential to determine the structure of genetic variation 
for a suite of traits within and across environments. 

Genetic variances depend on allele frequencies and are 
thus specific to populations and environments (Scheiner 
1993; Falconer and Mackay 1996). Heritability of a trait 
are typically affected by environmental conditions (e.g., 
Hoffmann and Merila 1999; Charmantier and Garant 
2005; Hallsson and Bjorklund 2012). This variation in 
heritability estimates across environmental conditions can 
be due to changes in the additive genetic variance (V A ) 
and/or the environmental variance (V E ), and may or may 
not entail a low genetic correlation between the expres- 
sions of the trait across the environments (Hoffmann and 
Merila 1999; Charmantier and Garant 2005). A priori pre- 
dictions about the influence of environmental conditions 
on the different components of heritability estimates are, 
however, hampered by discrepancies across studies (e.g., 
Hoffmann and Parsons 1991; Hoffmann and Merila 1999; 
Merila and Sheldon 1999; Charmantier and Garant 2005). 
The meta-analyses of studies comparing heritability esti- 
mates under favorable versus unfavorable conditions using 
data from wild populations by Charmantier and Garant 
(2005) suggested that, in general, estimates of both V A and 
V E are decreased under unfavorable conditions. The general 
trend in experimental research on Drosophila and some 
other insects under laboratory conditions, on the other 
hand, is that heritability estimates, including V A , are 
increased under more stressful conditions (e.g., Hoffmann 
and Merila 1999). It has been suggested that some of these 
differences may be because research under laboratory con- 
ditions often uses more extreme and/or more novel envi- 
ronmental stressors (Charmantier and Garant 2005). 
Importantly, the additive genetic variance has commonly 
been shown to increase under novel environmental condi- 
tions (independent of whether these conditions are favor- 
able or unfavorable), possibly due to expression of genes 
that have not been under selection in the more common 
environment (Holloway et al. 1990). Whereas organisms in 
the wild may typically experience stressful conditions, in 
the laboratory they do not. Therefore, once adapted to the 
optimal lab environment, introducing suboptimal condi- 
tions may in fact create novel environments. This highlights 
the importance of choosing a study system that facilitates 
disentangling these potentially confounding factors. One 
promising route might be to use a species for which it has 
been shown that adaptive phenotypic plasticity is an inte- 
gral part of its natural life history. Using such a study 
organism, the value and responses of genetic variance esti- 
mates of key traits in the predictable environment that 
drives the phenotypic plasticity and to an unpredictable 
stressful environment can be determined. Thus, within one 
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biological system, naturally perceived and more unpredict- 
able environments can be contrasted allowing the interpre- 
tation of the effects on genetic variation estimates in the 
light of past and future evolution. Here, we report on this 
approach using the tropical butterfly Bicyclus anynana. 

In B. anynana, phenotypic plasticity is a crucial com- 
ponent of the life cycle as it lives in highly seasonal envi- 
ronments for rainfall and temperature, and exhibits two 
very distinct seasonal forms that differ in wing pattern 
and many other traits (Brakefield et al. 2007). Environ- 
mental conditions during development, specifically those 
related to the thermal environment, are used as a cue for 
the future environment, and, subsequently, strongly influ- 
ence hormone dynamics, juvenile growth, and the result- 
ing adult life-history trajectories (e.g., Bauerfeind and 
Fischer 2005; Saastamoinen et al. 2010; Oostra et al. 
2011). Alterations in a suite of life-history and morpho- 
logical traits in B. anynana represent adaptive responses 
to seasonal differences in reproduction and survival as the 
wet and dry seasons are associated with favorable and 
more stressful environmental conditions, respectively (sea- 
sonal polyphenism; Brakefield and Larsen 1984; Brakefield 
et al. 2007, 2009). More specifically, individuals of the 
wet season form experience warmer ambient tempera- 
tures, and as a result have a shorter development time, 
become smaller as adults, reproduce at faster rate (higher 
investment to fecundity), and allocate less resources to 
body maintenance (i.e., fat reserves) resulting in reduced 
life span compared with the dry season form (e.g., Brake- 
field and Reitsma 1991; Brakefield and Kesbeke 1997; 
Pijpe et al. 2006). Dry season forms, on the other hand, 
tend to experience cooler ambient temperatures in the 
wild. Ambient temperature during the final larval instar is 
the main determinant of the two seasonal forms (Oostra 
et al. 2011). Recently, we have assessed in laboratory 
experiments how larval resources, which will also vary in 
nature, influence adult life-history traits in the wet season 
environment (Saastamoinen et al. 2010). Crucially in the 
context of our present study, even though developmental 
nutritional limitation generally reduces body mass and fit- 
ness (Bauerfeind and Fischer 2005; Saastamoinen et al. 
2010), individuals also changed their body allocation in 
ways likely to reflect an adaptive response to deteriorating 
environmental conditions (Van den Heuvel et al. 2013). 

Given that we have shown that phenotypic responses in 
B. anynana for a variety of traits work as adaptations allow- 
ing individuals to cope more effectively with one or other 
of the two alternating seasonal environments, responses to 
food limitation may also be season dependent. For 
instance, the increased thorax-to-abdomen ratio that 
occurs in females in response to larval food limitation in 
the wet season (Saastamoinen et al. 2010) may be less pro- 
nounced (or absent) in the dry season, as females in this 
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season already allocate very little to fecundity (i.e., resulting 
in higher thorax ratio; see Oostra et al. 201 1). It is therefore 
relevant to study plastic responses of the traits, and their 
interrelationships, under both wet and dry season environ- 
ments and to determine whether the potential for evolu- 
tionary change varies across seasons and environments. In 
particular the aim of this study was to test the following. 
First, as genotypes may be more constrained in reaching 
their potential under harsher environmental conditions 
(Gebhardt-Henrich and van Noordwijk 1991; Charmantier 
and Garant 2005), do we observe reduced heritability of 
performance traits in the laboratory-induced dry season 
form? Second, as nutritional limitation in the wild is less 
predictable and hence represents a more novel and less 
anticipated condition compared with the predictable ther- 
mal variation, do we observe increased heritability under 
conditions of nutritional limitation? Finally, as the levels of 
additive genetic variances can depend on their relative 
importance to fitness (Stearns and Kawecki 1994) and the 
strength or opportunity for selection, are the genetic vari- 
ance patterns consistent across the phenotypic traits? We 
implement an information-theoretic approach into a quan- 
titative analysis to arrive at an unbiased estimation of vari- 
ance components within and across environments that is 
insensitive to variation in trait means and measurement 
error. Thus, our study provides comprehensive estimates of 
genetic variation within and across environments in a spe- 
cies for which adaptive phenotypic plasticity is an integral 
part of its life history. We will discuss the results in the light 
of the effects of natural seasonal and unpredictable stressful 
environments and we will relate our findings to the ecology 
of Bicyclus. We conclude that although significant variation 
in genetic variance estimates exists between environments 
and traits, including for the strength of the across-environ- 
ment genetic correlation, no uniform pattern can be 
observed. The interpretation of heritability both for this 
study and in previously published literature will be greatly 
helped by studying variation in gene expression to allow a 
direct estimate of the absolute and relative contribution of 
these genes to the composite estimate that heritability is. 

Materials and Methods 
Study species 

Bicyclus anynana occurs in tropical and subtropical East 
Africa and feeds on fallen fruit as an adult. The butterflies 
used in this experiment originated from the stock popula- 
tion at the Leiden University, which was established in 
1988 from over 80 gravid females collected in Malawi. 
Several hundred butterflies are reared in each generation 
to sustain high levels of genetic variation (Van't Hof et al. 
2005). 



To ensure sufficient number of families for the experi- 
ment, over 70 mating pairs were established from the stock 
population. Two- to four-days-old virgin females were 
randomly mated with 2- to 6-day-old virgin males in an 
environmental chamber (+27°C, relative humidity 70%, 
L12:D12). Each female and male was allowed to mate only 
once. After mating, females were placed individually in a 
gauze-covered transparent pot with a young maize plant 
available for oviposition. Females were allowed to lay eggs 
for 5 days after which the eggs were removed and placed 
on a petri dish. Twenty-eight families with more than 80 
larvae were selected for the present experiment. 

Temperature and diet manipulation during 
development 

With the 28 families, we conducted a split brood experi- 
ment (Fig. 1) with two thermal environments (seasons), 
within which we further split each family into two larval 
food treatments. The first instar larvae in each family 
were split into two rearing temperatures; +27°C (RH 
70%, L12:D12) and +20°C (RH 70%, L12:D12) to mimic 
the natural environmental conditions during the wet and 
dry season, respectively (Brakefield 1997). Within each 
rearing temperature, each family was further divided into 
two groups of 20 larvae per plant kept inside a sleeve of 
gauze-like material to ensure nonstressful feeding densities 
for the family members. Sleeves were checked daily and 
fresh plants were provided when needed to ensure ad libi- 
tum feeding for the larvae. 

On the first day of the final, fifth, instar, the larvae within 
each thermal treatment were further randomly assigned 
into two larval food treatments: food limitation and con- 
trol. Equal numbers of individuals within each family were 
assigned to each treatment. Larvae were transferred individ- 
ually to a petri dish with either fresh maize leaves (control) 
or a piece of set agar (1.5 cm 2 to ensure humidity; food 
limitation). The larvae were kept in the petri dishes (with 




Figure 1. Experimental design: Twenty-eight full-sib families (A) were 
each split and reared under both dry and wet season conditions (6; +20 
and +27°C, respectively). In the final instar, family groups from both 
thermal conditions were further split into control and food limitation 
treatment (Q. As individuals eclosed six females and six males from 
each separate family group were assessed for the five life-history traits (D). 
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or without food) for 2 or 3 days in the warm (wet season) 
and cooler (dry season) conditions, respectively. One extra 
day of food limitation was added under the cooler thermal 
conditions in order to have equal reduction in the body 
mass due to the food limitation (based on a pilot experi- 
ment; M. Saastamoinen, unpubl. data), as larvae, in general, 
develop more slowly under cooler temperatures. After the 2 
or 3 days on a petri dish, individuals were transferred back 
onto potted maize plants (10 individuals per plant and fami- 
lies separate) with subsequent ad libitum feeding. Larvae 
from the two different treatment groups and from the differ- 
ent daily cohorts were kept separate to enable development 
time to be recorded for each individual. Individuals were 
sexed and weighed 1 day after pupation (Sartorius micro- 
balance GMBH; Sartorius AG, Goettingen, Germany), after 
which individuals were placed individually in small con- 
tainers in which the resulting adult butterflies eclosed. 

Life-history traits measured 

We assessed development time, pupal mass, resting meta- 
bolic rate (RMR), thorax ratio, and fat percentage for three 
males and three females from each sleeve (i.e., whenever 
possible six females and six males/family/treatment/sea- 
son). RMR was measured 1 day after eclosion with a Li-Cor 
LI-6251 C0 2 analyzer in a Sable systems respirometer setup 
with a push-through flow of 100 mL/min. RMR of the 
individual butterflies was measured in small cylindrical 
glass containers, which were kept in a temperature-con- 
trolled climate chamber (+27°C and +20°C, respectively, 
for wet and dry season butterflies) for three consecutive 
measurements of 12 min each. Butterflies were measured 
in the dark to ensure a resting state of the butterfly. The 
first C0 2 reading of each individual was discarded and the 
second and third readings were averaged before analysis to 
only include measurements from inactive butterflies. C0 2 
data from the two remaining consecutive measurements 
were then analyzed using Datacan 5.4 (Sable system, Berlin, 
Germany). Butterflies were sacrificed after measurement of 
RMR and the abdomen and thoraces were dried (60°C for 
24 h) and weighed (Sartorius microbalance GMBH; Sarto- 
rius AG) to assess thorax ratio (thorax dry mass/(thorax 
dry mass + abdomen dry mass)). For the fat percentage 
analyses, the dried thoraces and abdomens were submerged 
in ethyl acetate for 4 days, dried at 60°C for 24 h, and 
weighed again (Brakefield et al. 2009). The difference 



between initial dry weight and weight after fat extraction 
was the absolute fat content. Fat percentage was calculated 
as the ratio of the absolute fat content to the total dry mass. 

Data analyses 

Assessment of phenotypic differences in the 
response to season and food limitation 

A linear mixed model approach (SAS version 9.2 for Win- 
dows; SAS Institute Inc., Cary, NC) was used to examine 
the effects of temperature and food treatments on pheno- 
typic variation. In all the analyses, sex was included as a 
fixed factor, and family and the sleeve nested within a fam- 
ily were included as a random factor. We used backward 
model selection by starting with a full model for each trait 
and sequentially eliminating terms with lowest f-values 
until all terms in the model were significant. Development 
time, fat percentage, and thorax ratio were log transformed 
so that values were normally distributed. In the analyses of 
the RMR fat-free dry mass was used as a covariate. 

Assessment of genotypic differences in the 
response to season and food limitation 

All families were reared under four different environmen- 
tal conditions mimicking the wet or dry season (tempera- 
ture), and with or without food limitation. In the 
analyses, the traits measured were viewed as specific to 
each of these four environments (character state; Lynch 
and Walsh 1998). Because we used a full-sib design, we 
could only estimate genetic (co)variances, which include 
- apart from additive genetic effects - dominance and 
maternal effects. Thus, the heritability we estimated must 
be considered as the upper estimate under the assump- 
tions that dominance and maternal variance components 
are zero. Family members were reared in different com- 
mon environments (sleeves), thereby minimizing resem- 
blance across relatives due to rearing effects. Each trait 
has an environment-specific genetic variance with genetic 
covariances between the environment-specific trait expres- 
sions. Thus, there are in total four genetic variance {V G ) 
and six genetic covariances (Cg) between all combinations 
of food limitation (n = no stress; s = stress) and season 
(wet and dry defined by thermal conditions). The (co) 
variances can be summarized as the matrix, 



V G (n,wet) 
C G («, wet — s, wet) V G (s, wet) 

C G (n, wet - n, dry) C c (s, wet - n, dry) V G (n, dry) 

_ C G (n, wet - s, dry) C G (s, wet - s, dry) C G (n, dry - s, dry) V G (*;,dry)_ 
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The genetic covariances can be scaled to a genetic cor- 
relation following the standard definition, for example, 
r G (n,wet — s, wet) = C G (n,wet — s, wet)H{V G {n, 
wet) x V G (s,wet)). Because we estimate compounded 
genetic (co)variances (see above), this matrix is an 
approximation of the G matrix with additive genetic (co) 
variances. Whereas the V G estimate can be considered as 
an upper estimate for the additive genetic variance, the 
C G estimate lacks mathematical relationship to the addi- 
tive genetic covariances because the covariances on the 
numerous levels (additive genetic, dominance, maternal 
effect) compounded in C G need to align. Hence, our esti- 
mates of Cg and thus of r G must be interpret with cau- 
tion and are not therefore the main focus of this study. 

Interactions between genotype and environment (GEI, 
defined here by the unique combination of season and 
food limitation) can take two nonmutually exclusive 
forms (Lynch and Walsh 1998). First, the relative ranking 
of the breeding values may change between environments. 
In general, a negative or low genetic covariance between 
environments indicates that the ranking of genotypes is 
changed (crossing reaction norms). Second, additive 
genetic variances may be specific to the environment. 
However, changes in additive genetic variances across 
environment are likely to be subjected to scaling, where 
not only the additive genetic variance changes but also 
the residual variance (and hence the phenotypic variance). 
Approaches to study changes in genetic architecture inde- 
pendently from scaling include standardization of vari- 
ance components with the traits mean (coefficient of 
variation) or standardization of trait values with their var- 
iance (variance standardization; Lynch & Walsh 1998). 
Here, we follow the latter approach and thereby explore a 
form of GEI which creates changes in heritability across 
the four environments. The G matrix (eq. 1) was esti- 
mated by defining the linear mixed model, 

y = X/J + Zu + e, (2) 

where y is a vector of observations on all individuals, ft is 
a vector of fixed effects, X represents a design matrix (of 
Os and Is) relating to the appropriate fixed effects to each 
individual, u is a vector of random effects, Z is a design 
matrix relating the appropriate random effects to each 
individual, and e is a vector of residual errors. G is 
defined as the matrix for vector u, and its elements 
(the genetic (co)variances) can be estimated by using 
information on the coefficient of coancestry 0 { j between 
individuals i and j, which is directly obtained from the 
pedigree. All individuals measured were the descendants 
of butterflies mated in a full-sib cross. There were 56 base 
parents with a total of 1206 descendants. The genetic 
effects in environment E (i.e., the combination of food 
limitation or no food limitation, and dry or wet season 



conditions experienced) were assumed to be normally dis- 
tributed with mean of zero (i.e., defined relative to the 
environment-specific fixed-effect mean) and with an 
genetic variance of (T 2 ^. This variance (and the additive 
genetic covariance between all E) was estimated by REML 
from the variance-covariance matrix of additive genetic 
effects which is equal to Act 2 , where A has elements, 

Ajj = 2&ij. (3) 

The fixed-effect structure of equation (2) accounted for 
variation in age of the butterflies when they entered the 
final instar and for variation between the sexes. We thus 
considered all data on both sexes in order to maximize our 
power to detect changes in genetic variances across envi- 
ronments. By fitting "sex" as a fixed effect, we corrected 
only for the difference in the mean trait expression between 
the sexes and thereby assumed that the between-sex genetic 
correlations for traits did not differ from +1 (no gene- 
by-sex interaction). In addition, for RMR we included the 
total fat-free dry weight of the individual as a fixed effect. 
Residuals were assumed to be heterogeneous (environment 
specific) and not correlated across environments. 

Variances in a linear mixed model are conditional 
upon the fixed-effect structure. Mixed model phenotypic 
variance is, in this case, the sum of the REML genetic 
(including dominance and maternal variances) and resid- 
ual variances. We incorporated variance scaling by stan- 
dardizing the raw data prior to analysis to have a REML 
variance of unity (1) in each environment. This was done 
by first running a model that only included the fixed 
effects, where the four residual variances (assumed to be 
uncorrelated across environments) estimated the environ- 
ment-specific REML variances. In further analyses, the 
data were divided by the environment- specific REML 
standard deviation. By doing so, all trait values become 
dimensionless (expressed in unit REML phenotypic SD) 
and the diagonal in equation (1) thereby consisted of the 
upper estimates of the trait-specific heritabilities. 

Given the four environments considered, there are, for 
each trait, 15 models to consider and, in addition, the null 
model with residuals only (no heritability). Models were 
implemented in ASReml (VSN International, Hemel 
Hempstead, UK), which provides the log likelihood of the 
mixed model. Because many of the models are not nested, 
model comparison relied on an information theoretical 
approach based on the Akaike information criterion (AIC; 
Akaike 1974; Wagner et al. 1997; Burnham and Anderson 
2002). AIC was calculated as -2log(L) + 2K, where log(L) 
was the model's log likelihood and K the number of param- 
eters estimated. All models were ranked in ascending order 
based on their AIC, where a difference in AIC of more than 
two compared to the model with the lowest AIC was con- 
sidered as evidence of deterioration in model fit (Burnham 
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and Anderson 2002). In ASReml, constraining the diagonal 
of the G matrix (eq. 1) does not constrain the genetic co- 
variances, which were left unconstrained in all models. We 
therefore did not consider the covariances or the fixed 
effects in calculating K because these parameters were esti- 
mated in all models and hence are factored out when doing 
model comparisons based on AIC. Thus, we calculated K as 
the number of genetic variances estimates, and K ranged 
from 0 (model with residuals only) to 4 (environmental- 
specific variances). Akaike weights w for model i was calcu- 
lated as Wi = exp(AAIQ)/£exp(AAIC), where AAIQ is the 
difference in AIC between model i and the top model (i.e., 
the model with the lowest AIC). Models ranked within two 
AIC units of the top model were considered as reasonable 
candidate models (Burnham and Anderson 2002). Because 
one typically finds some level of support for multiple candi- 
date models, model averaging is advocated to provide more 
precise estimates (Burnham and Anderson 2002). We 
model averaged the estimates of the genetic and residual 
variances across all 15 models, where the model-averaged 
variance V* was calculated by weighing V b the variance 
estimate of model such that V* = E(Vi x wi) with its 
model-averaged standard error SE* calculated as 
SE* = ]T[w; x ^(SE 2 + (Vi - V*) 2 )] (Burnham and 
Anderson 2002). 

Results 

Phenotypic variation 

Larvae developed faster and were smaller as pupae under 
the warmer wet season conditions compared to the cooler 
dry season conditions, and males developed faster and 
were smaller as pupae than females (Table 1, Fig. 2A and 
B). For both of these traits, sex differences were slightly 
larger in dry season conditions (sex*season interaction; 
Table 1). Experiencing food limitation during the last 
instar increased development time and decreased pupal 
size in both seasons (Table 1, Fig. 2A and B). The effect 
on development time was larger in the dry season (food 
limitation*season interaction; Table 1): food-limited larvae 
took, on average, an extra 4 and 6 days to pupate in wet 
and dry season conditions, respectively. Importantly, 
however, the length of the food treatment (i.e., days with- 
out food) was also longer in the dry season due to the 
generally longer development time under these condi- 
tions, and hence in both seasons the food-limited individ- 
uals needed twice the number of days they had been 
without food in the experiment to complete development 
into an adult butterfly (for more details see Material and 
Methods). The effect of food limitation on pupal mass on 
the other hand was significantly stronger in the wet 
season (food limitation*season interaction; Table 1 and 



Fig. 2B), as the decrease was 4% and 7%, respectively, for 
dry and wet season. 

Thorax ratio was higher in males than in females 
(Table 1, Fig. 2C). The effect of season on thorax ratio 
was significant (Table 1), however, a significant interac- 
tion between sex and season indicated that the relative 
thorax ratio was only higher in females under dry season 
conditions (Table 1 and Fig. 2C; see also Oostra et al. 
2011). Experiencing stressful food-limited conditions dur- 
ing larval development increased relative allocation to the 
thorax in both sexes and in both seasons (Table 1, 
Fig. 2C). The change in the relative allocation to the tho- 
rax is caused by a decrease in both abdomen and thorax 
mass due to food limitation, with a more marked 
decrease in the former (for females: 10% and 8%; and for 
males: 11% and 6% decrease, respectively, for abdomen 
and thorax mass in relation to food limitation). 

Males had a higher fat percentage than females (Table 1, 
Fig. 2D). For males the fat percentage was higher in the wet 
season compared with the dry season, whereas for females 
the opposite occurred, and a higher fat percentage was 
found in this sex under dry season conditions (sex*season 



Table 1. Analyses of the effects of season, feeding treatment, and 
sex on phenotypic traits measured. In all the analyses, family and 
sleeve nested within family were included as random factors. 





df 


F 


P 


Development time 








Season 


1,1084 


12642.2 


<0.0001 


Feeding treatment 


1,1084 


6440.1 


<0.0001 


Sex 


1,1084 


215.2 


<0.0001 


Season*Sex 


1,1084 


9.4 


0.0022 


Season*Feeding treatment 


1,1084 


124.5 


O.0001 


Pupal mass 








Season 


1,1079 


197.4 


<0.0001 


Feeding treatment 


1,1079 


104.0 


O.0001 


Sex 


1,1079 


2795.4 


<0.0001 


Season*Sex 


1,1079 


7.6 


0.006 


Season*Feeding treatment 


1,1079 


9.1 


0.003 


Thorax ratio 








Season 


1,1090 


15.5 


O.0001 


Feeding treatment 


1,1090 


31.5 


<0.0001 


Sex 


1,1090 


10361.9 


<0.0001 


Season*Sex 


1,1090 


30.8 


<0.0001 


Fat percentage 








Season 


1,1087 


1.7 


0.187 


Feeding treatment 


1,1087 


13.7 


0.0002 


Sex 


1,1087 


1371.8 


<0.0001 


Season*Sex 


1,1087 


19.8 


<0.0001 


Season*Feeding treatment 


1,1087 


18.8 


<0.0001 


Resting metabolic rate 








Season 


1.10S8 


301.9 


<0.0001 


Feeding treatment 


1,1058 


4.8 


0.029 


Sex 


1,1058 


55.9 


<0.0001 


Fat-free mass 


1.10S8 


139.5 


O.0001 
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interaction; Table 1 and Fig. 2D). Experiencing nutritional 
limitation decreased fat percentage in general (Table 1). 
However, there was a significant interaction between the 
season and effect of stress: food limitation had a negative 
effect on fat percentage in the dry season, but no effect in 
the wet season (Table 1 and Fig. 2D). 

Heavier individuals had a higher RMR, yet after taking 
into account the dry body mass, males had a higher RMR 
than females (Table 1). RMR was additionally higher in 
individuals reared under wet season conditions (Table 1), 
but this difference is largely due to the ambient tempera- 
tures during the measurements, which were warmer in 
wet than in the dry season conditions (see Material and 
Methods). Individuals that had experienced nutritional 
limitation had a lower RMR than when fed ad libitum 
throughout their development (Table 1 and Fig. 3E). 

Genetic variation and genotype- 
environment interactions 

The difference in AICs when comparing the model with 
one constant heritability (model "AAAA") to the model 
with residuals only was at least 29 AIC units (Table 2), pro- 
viding clear evidence that heritable variation was present 



Figure 2. Effects of season (solid and dashed 
line for dry and wet season, respectively) and 
larval feeding treatment on development time 
(A), pupal mass (6), thorax ratio (0, fat 
percentage (D), and resting metabolic rate 
(RMR; corrected for fat-free dry mass), (£) 
plotted separately for females (red line) and 
males (blue line). Error bars represent SE of 
means. 

for all traits in this population. Differences in heritability 
across the different environments were, however, less clear. 
Typically, there were several candidate models within two 
AIC units of the top-ranked model (Table 2). We thus 
investigated these top models for consistent patterns. For 
pupal mass and development time, there was no strong sta- 
tistical evidence of any GEI, as models with equal heritabili- 
ties in all environments (model "AAAA") were within two 
AIC from the best model (Table 2, Fig. 3). 

For thorax ratio, we found good evidence that the dry 
season increased heritability independently of whether the 
developing larvae were put under food limitation or not 
(Table 2, Fig. 3). The covariances between all these envi- 
ronments were strongly positive (Table 3) indicating that 
ranking of the genotypes was generally maintained. 
Hence, we find a heritability-increasing genotype - dry 
season interaction for thorax ratio. 

The generally low heritability in RMR makes inferences 
difficult. We find no strong evidence for GEI, as the 
model with equal heritability in all four environments 
(model "AAAA") falls within the top models (Table 2). 
In addition, the heritability under food limitation in the 
dry season tended to be negative (Fig. 3, Table SI), indi- 
cating difficulties in estimating the environment-specific 
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Figure 3. Model-averaged restricted maximum 
likelihood (REML) estimates of environment- 
specific upper estimate of heritability with its 
standard error for development time (A), pupal 
mass (B), thorax ratio (0, fat percentage (D), 
and resting metabolic rate (RMR, E). Four 
environments are characterized by food stress 
(standard vs. food limitation) and season (wet 
season vs. dry season) the larvae experienced 
during development. Because the data were 
standardized prior to analysis to have REML 
variance of unity, the plotted values 
approximate the environment-specific 
heritabilities (i.e., V G + V R rj 1 in each 
environment). Model averaging was performed 
across all 15 candidate models shown in 
Table 2. For the traits "RMR" and "fat 
percentage", the REML genetic variances were 
negative in the dry season under food 
limitation and in the wet season under food 
limitation, respectively. Further details on the 
models' fixed effects, residual variances, and 
estimates are provided in the Supplement. 
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heritability. The largest heritability is found in the dry 
season without food limitation (Fig. 3), but the estimate 
is still very modest (0.094). 

The pattern of heritability in fat percentage showed 
similarities to that of thorax ratio in that there was a clear 
increase in heritability in the dry season when larvae 
developed under optimal food conditions (Table 2 and 
Fig. 3). In addition, there was some evidence of geno- 
type-food limitation interactions for fat percentage in 
both seasons as many of the top models included this 
interaction (Table 2, Fig. 3). However, for the wet season 
this interaction was mainly driven by the estimate of heri- 
tability turning negative in the "food limitation/wet" 
environment (see Table S2 for alternative models where 
this component is constrained), and this genotype-food 
limitation interaction was thus not well supported. Fat 
percentage was the only trait where the genotype-food 
treatment-season interaction model (model "ABCD" with 
specific heritability in each environment) was ranked 
within two AIC from the top model, indicating substan- 
tial plasticity in this trait. Although interpretation is 



hampered by the apparent absence of genetic variance in 
the wet season under food limitation, the genetic relation- 
ships between the environments with non-negative herita- 
bility were clearly positive (Table 3), showing that ranking 
of genotypes for this trait across the environments was 
generally maintained. 

Discussion 

Changes in environmental conditions are known to influ- 
ence quantitative genetic variation, but there are some 
discrepancies among studies (e.g., Hoffmann and Merila 
1999; Charmantier and Garant 2005). In this study we 
specifically aimed to assess phenotypic and genetic 
(genetic variances and genetic covariances) responses to 
commonly experienced and predictable environmental 
fluctuation, and to compare them with responses to a less 
predictable environmental fluctuation. Thermal condi- 
tions during development represented the former 
conditions as it is an environmental variable that deter- 
mines the seasonal polyphenism in the study organism 
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Table 2. Comparison of models allowing for different additive genetic variance structures in the G matrix. 
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Model 


log(Z.) 


AAIC 


Weight 


Model 


logtt) 


AAIC 


Weigl 


npvplnnmpnt timp 








Pi in^l 

r uugi ii i ci _> _> 








AABA 


— 590.00 


0.00 


0.23 


AAAB 


—546.42 


0.00 


0.23 


AAAA 


— 591 .10 


0.20 


0.18 


AAAA 


—547.49 


0.13 


0.20 


AAAB 


— 590.23 


0.47 


0.14 


ABAB 


—546.52 


0.20 


0.19 


ABAB 


— 590.40 


0.80 


0.10 


AABA 


—547.1 1 


1 .36 


0.06 


ABBA 


— 590.51 


1 .02 


0.08 


ABBA 


—547.22 


1 .60 


0.05 


AABC 


— 589.66 


1 .34 


0.06 


ABBB 


—547.25 


1 .66 


0.04 


ABCB 


— 589.86 


1 .73 


0.04 


ABAC 


—546.29 


1 .73 


0.04 


ABCA 


-589.97 


1 .94 


0.03 


AABB 


-547.32 


1 .79 


0.04 


AABB 


— 591 .10 


2.20 


0.03 


ABBC 


-546.38 


1 .92 


0.03 


ABBB 


— 591 10 


2 20 


0 03 


AABC 


-546.41 


1 97 


0 03 


ABAA 


— 591 10 


2 20 


0 02 


ABAA 


-547.49 


2 1 3 


0 03 


ABAC 


— 590 1 0 


2 21 


0 02 


ABCB 


-546.52 


2 20 


0 03 


ABBC 


— 590 1 4 


2 30 


0 02 


ABCA 


-547.1 


3 34 


0 01 


ABCD 


— 589 65 


3 31 


0 01 


ABCC 


-547.24 


3 62 


0 01 


ABCC 


— 591 10 


4 20 


0 00 


ABCD 


-546.29 


3 72 


0 01 




—626 34 


68 76 




res 


-626.08 


1 55 49 




Thnrav ratin 








Resting metabolic rate 






AABB 


— 593.90 


0.00 


0.53 


AABA 


-616.14 


0.00 


0.31 


AABC 


— 593.78 


1 .76 


0.09 


ABAB 


-616.42 


0.55 


0.18 


ABCC 


— 593.89 


1 .98 


0.07 


AAAA 


-617.76 


1 .23 


0.09 


AAAA 


— 595.91 


2.03 


0.07 


AAAB 


-616.79 


1 .29 


0.08 


AABA 


— 594.99 


2.18 


0.06 


AABC 


-615.89 


1 .49 


0.07 


ABBB 


— 595.05 


2.30 


0.05 


AtSLb 


(Tir QT 

— d i b.y j 


1 .57 


0.06 


ABAA 


— 595.36 


2.93 


0.03 


ABBA 


-617.01 


1 .73 


0.06 


ABCB 


-594.38 


2.96 


0.03 


ABCA 


-616.14 


2.00 


0.04 


ABCD 


-593.76 


3.71 


0.01 


ABAC 


-616.29 


2.30 


0.03 


AAAB 


— 595.77 


3.74 


0.01 


AABB 


-617.58 


2.88 


0.02 


ABBA 


— 595.78 


3.76 


0.01 


ABAA 


-617.59 


2.88 


0.02 


ABCA 


— 594.82 


3.84 


0.01 


ABBC 


-616.67 


3.05 


0.01 


ABAB 


— 595.89 


3.98 


0.01 


ABBB 


-617.75 


3.22 


0.01 


ABBC 


— 595.05 


4.30 


0.01 


ABCD 


-615.83 


3.37 


0.01 


ABAC 


— 595.28 


4.77 


0.00 


ABCC 


-617.54 


4.79 


0.00 




—626.30 


61 .42 




res 


-610.60 


30.65 




pat nprrpnt^np 
















ABCA 


— 599 25 


0 00 


0 51 










ABCC 


— 599 73 


0 97 


0 1 9 










ABCD 


— 598 95 


1 40 


0 13 










AABB 


—601 30 


2 09 


0 06 










ABAA 


—601 32 


2 1 5 


0 06 










AABC 


600 66 


2 82 


0 03 










ABAC 


—601 21 


3 93 


0 01 










AABA 


-602.29 


4.08 


0.01 










ABCB 


-602.27 


6.04 


0.00 










ABAB 


-604.46 


8.41 


0.00 










AAAA 


-606.18 


9.87 


0.00 










ABBB 


-605.50 


10.50 


0.00 










ABBA 


-605.96 


11.42 


0.00 










AAAB 


-606.13 


11.76 


0.00 










ABBC 


-605.50 


12.50 


0.00 










res 


-625.29 


44.75 













Model structure follows the diagonal of equation (1) and letter coding refers to the (dis)similarity of the additive genetic variance for (left to right) 
no food limitation in the wet season, food limitation in the wet season, no food limitation in the dry season, and food limitation in the dry sea- 
son. Same letters indicate that the respective additive genetic variances were constrained to be the same. The model with only residuals is indi- 
cated as "res". For each model, the log likelihood (log(A)) is presented. All models are ranked in ascending order for their AIC value and, for each 
model, the difference in AIC value with the model with the lowest AIC value is presented (AAIC). All models with an AIC value within 2 units 
from the model with the lowest AIC value are printed in bold and are all considered as reasonable descriptions of the data. 
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Table 3. Model-averaged genetic covariances with their standard error and (in brackets) the scaled covariance (genetic correlation) based on the 
model-averaged covariances and variances (variances reported in Table S1 and Fig. 3). 

Wet season/No food limitation Wet season/Food limitation Dry season/No food limitation 

Development time 



Wet season/Food limitation 


0.23 ± 0.090 (0.86) 






Dry season/No food limitation 


0.13 ± 0.080 (0.56) 


0.0042 ± 0.082 (0.018) 




Dry season/Food limitation 


0.13 ± 0.090 (0.44) 


0.15 ± 0.092 (0.50) 


0.10 ± 0.083 (0.50) 


Pupal mass 








Wet season/Food limitation 


0.35 ± 0.12 (0.88) 






Dry season/No food limitation 


0.37 ± 0.12 (0.99) 


0.34 ±0.12 (0.88) 




Dry season/Food limitation 


0.39 ± 0.13 (0.91) 


0.40 ±0.14 (0.88) 


0.44 ±0.13 (1.0) 


Thorax ratio 








Wet season/Food limitation 


0.08 ± 0.06 (0.71) 






Dry season/No food limitation 


0.08 ± 0.07 (0.04) 


0.15 ± 0.07 (0.77) 




Dry season/Food limitation 


0.15 ± 0.07 (0.86) 


0.09 ± 0.07 (0.49) 


0.26 ±0.10 (0.92) 


RMR 








Wet season/Food limitation 


0.02 ± 0.04 (1.1) 






Dry season/No food limitation 


0.07 ± 0.05 (1.4) 


0.004 ± 0.05 (0.13) 




Dry season/Food limitation 


N/A 


N/A 


N/A 


Fat percentage 








Wet season/Food limitation 


N/A 






Dry season/No food limitation 


0.18 ± 0.08 (0.92) 


N/A 




Dry season/Food limitation 


0.11 ± 0.07 (0.79) 


N/A 


0.14 ± 0.08 (0.56) 



Covariances and correlations are presented in matrix form between all four environments. Environments are wet season or dry season/food limita- 
tion or No food limitation. Note that the variance-covariance matrix was not constrained to be positive definite and some of the scaled covari- 
ances therefore fall outside the interval of -1 to 1. For the traits "resting metabolic rate (RMR)" and "fat percentage", the REML genetic 
variances was negative in the dry season under food limitation and in the wet season under food limitation, respectively (Fig. 3, Table S1), and 
hence genetic covariances could thus not be defined (reported here as "N/A"). 



B. anynana (i.e., wet vs. dry season morph which differs 
in both wing morphology [data not shown] and life his- 
tory; e.g., de Jong et al. 2010; Oostra et al. 2011) and 
hence for which the species has a past history of selection. 
Food availability during development, on the other hand, 
was chosen to represent the nonpredictable and hence 
more novel environmental variable. 

Season-dependent responses to food 
limitation 

As expected, the thermal conditions experienced during 
larval development produced the two distinct seasonal 
forms. Hence, individuals under the wet season thermal 
environment developed faster, were smaller, had a smal- 
ler thorax ratio, and lower fat percentage (the last two 
for females only). The food limitation during develop- 
ment had an equal effect between the two seasons for 
some traits, whereas for others the responses were sea- 
son specific. For development time, thorax ratio, and 
RMR, the responses to food limitation were equal 
between wet and dry seasonal environments. The signifi- 
cant interaction between season and food limitation for 
development time is explained by the initial difference 
of food-limited days (see Material and Methods), as in 



both cases the development time of the food-limited 
larvae increased by about twice the number of food- 
restricted days they experienced as larvae. For thorax 
ratio, the nonseason-specific response to food limitation 
is somewhat in contradiction to our recent theoretical 
model, which showed how an increase in thorax ratio 
due to poor developmental conditions could be adaptive 
in the wet season, as only in this season high-quality 
habitats exist and hence can be located by individuals 
via dispersal (Van den Heuvel et al. 2013). Thorax ratio 
can, however, conversely also measure allocation to 
reproduction, as a low thorax ratio indicates that indi- 
viduals have allocated more resources into the abdomen 
(i.e., eggs in females) at the cost of a relatively small 
thorax. Therefore, the smaller thorax ratios under more 
optimal conditions may indicate increased investment in 
reproduction under these conditions, whereas a higher 
thorax ratio in the dry season may reflect increased 
allocation to dispersal/food searching capacity. Finally, 
alteration of metabolism occurs in other organisms (e.g., 
Caenorhabditis elegans; Van Voorhies 2002) in response 
to rapid environmental change including an absence of 
food, and presumably decreased RMR under food limi- 
tation can enable individuals to cope better under 
stressful conditions. 



2013 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



3585 



Quantitative Genetics of Phenotypic Plasticity 

The effect of food limitation on pupal mass was larger 
under the wet season environment. This may indicate that 
the larvae under dry season conditions are better adapted 
to compensate for food limitation. This is reasonable as 
they are more likely to experience such conditions in the 
wild. Alternatively, nonadaptive explanations may relate 
to constraints of the thermal environment on physical 
processes of growth and resources acquisition and alloca- 
tion (e.g., van der Have and de Jong 1996). The effect of 
food limitation on fat percentage was similarly season 
dependent, as food limitation decreased fat percentage 
only in the dry season. The absolute fat content is higher 
in dry season butterflies as a result of their larger size, 
and the relative fat content has been shown to be higher 
as well over a range of temperatures spanning from the 
dry to the wet season condition (Oostra et al. 2011; but 
see our Fig. 2). Fat in the dry season is likely to be used 
for maintenance and needed for increased stress resistance 
and increased life span in that season (Pijpe et al. 2007). 
Our interpretation is that allocation to fat reserves, and 
hence surviving, under dry season condition is probably 
maximized, and reduced resources during the larval stage 
will therefore greatly impact fat content. In contrast, lar- 
vae under the wet season conditions are likely to allocate 
fat to reproduction, and the increased fat allocation under 
food limitation may be a predictive response to the forth- 
coming stressful conditions potentially allowing individu- 
als to cope with harsher future conditions (Saastamoinen 
et al. 2010). 

Changes in heritability across 
environmental conditions 

We used a full-sib design to estimate the upper value of 
heritability. Our estimate is the upper value because our 
design cannot partition out genetic dominance variance 
and maternal effects from the additive genetic variance. 
We have standardized the variances prior to analyses in 
order to be sure that changes in the genetic architecture 
are not merely due to scaling. Variance standardization 
means that we focus on whether heritability estimates of 
the traits change across environments. We feel that 
changes in heritability provide not only an intuitive level 
to study GEI, but heritabilities also provide the best pre- 
diction for the potential of evolutionary change caused by 
natural selection. The drawback of variance standardiza- 
tion is that an equal heritability across two or more envi- 
ronments does not equate to equal genetic variances, but 
rather states that differences in variances between envi- 
ronments affect both genetic and environmental variances 
to the same degree (Houle et al. 2011). In addition, 
because many non-nested models are compared, we have 
used model averaging to obtain precise estimates that take 
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into account model uncertainty (Burnham and Anderson 
2002). 

We find significant heritability for all the life-history 
traits, although for RMR the estimate was low. In general, 
the traits we considered appear to show little GEI: only in 
two of the five traits (thorax ratio and fat percentage) 
heritabilities changed across thermal environment. For fat 
percentage there was additionally an indication of GEI 
across the nutritional environment. The season-specific 
heritabilities were not in line with our predictions. We 
predicted, following Charmantier and Garant (2005), 
reduced heritability in the dry season environment (low 
ambient temperature), as genotypes may be more con- 
strained under harsher environmental conditions. How- 
ever, for both thorax ratio and fat percentage, the 
heritability increased under a dry season environment. 
Furthermore, we assumed that food limitation would rep- 
resent an unpredictable environment in our study system, 
and hence predicted increased heritability under such 
conditions, especially in the wet season (e.g., Messina and 
Fry 2003; Dmitriew et al. 2010). Genetic variance did not, 
in general, react to this change in environmental condi- 
tions, with the exception of fat percentage, where - con- 
trary to our prediction - heritability tended to decrease 
under food-limited conditions. 

Apart from causing changes in the additive genetic var- 
iance, GEI can also operate by uncoupling the genetic 
correlations between environments (Falconer and Mackay 
1996). Generally, for most of the traits we estimated posi- 
tive genetic covariances, although our estimates should be 
interpreted with some caution as our design was not opti- 
mal for assessing these (see Material and Methods). Nev- 
ertheless, our findings suggest that the main source of 
GEI in the B. anynana system concerns differential 
expression of genetic variances, rather than low or nega- 
tive genetic correlations between environments. Hence, it 
seems likely that the phenotypic values for these traits are 
governed by the expression and regulation of the same 
genes and pathways in the different environments (Falconer 
and Mackay 1996). However, the reported correlations 
leave room for season-specific variation and this may in 
part be related to the number of genes contributing to the 
trait in each season which may explain the variation in 
heritability between the seasons. 

As the levels of genetic variance in a trait are predicted 
to depend on the relative importance of the trait to fit- 
ness (Stearns and Kawecki 1994) and the strength/oppor- 
tunity for selection, we were interested in whether the 
heritability patterns were consistent across the phenotypic 
traits. As stated above, similar season-specific patterns 
were estimated between thorax ratio and fat percentage. 
Additionally, for development time and pupal mass, we 
estimated relatively high heritability, but little variation 
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across environments. Our findings thus underline that 
generalization of patterns of GEI (Charmantier and 
Garant 2005) is difficult, especially when dealing with 
different types of environmental condition (natural vs. 
experimentally imposed). 

Interpretation of the results in the light of 
the evolutionary history and ecology of 
Bicyclus 

First, we need to address that not finding statistical evi- 
dence for an interaction (GEI in our case) may be due to 
lack of power to detect this interaction. Unfortunately, 
the power of our quantitative genetic approach is not eas- 
ily calculated. It should be noted, however, that our 
design had sufficient power to detect overall significant 
levels of genetic variances for all traits, including variable 
levels per trait, from low for RMR to high for pupal 
mass. This is generally in line with quantitative genetic 
and life-history theory (Stearns 1992; Roff 1997). The 
model-averaged estimates of heritability for developmen- 
tal time and pupal mass in each food/seasonal environ- 
ment are strikingly close, suggesting that constancy of 
heritability across the different environments is the most 
likely cause of not finding GEI in these traits. For RMR, 
on the other hand, there simply is very little overall 
genetic variance and thus little "room" for interactions 
with environment. 

Our finding that heritability is higher in the dry season 
for thorax ratio and fat percentage could potentially be 
explained by differential selection pressures between the 
seasons. Reduced larval food sources at the transition 
between the wet and the dry season and unpredictable 
resources for the adult in the dry season is part of the 
environmental variation the species naturally experiences. 
As thorax ratio and fat percentage are part of the adaptive 
suit of traits, these traits have been subjected to natural 
selection. We would therefore expect lower heritability for 
these traits in a dry season environment. However, we 
have measured variation in these traits at eclosion, which 
may not be the time in the life history when these traits 
contribute the most to fitness, and thus the time when 
the traits are under the strongest selection. Moreover, 
even though fat percentage in the dry season is strongly 
linked with adult starvation resistance in that season, 
selection on thorax ratio and fat percentage may also 
operate in the wet season, which is the main reproductive 
season of the butterfly. As the heritability in a given trait 
is predicted to depend on its relative contribution to fit- 
ness (Stearns and Kawecki 1994) and the strength/oppor- 
tunity for selection, this would be according to the 
general predictions. However, as we lack measures of 
selection operating on life-history traits in the various 



seasons, in contrast to estimates for wing pattern varia- 
tion (Brakefield and Frankino 2009), such a scenario 
remains speculative. Moreover, selection is likely to be 
complex in the dry season as it must favor an effective 
and well-timed switch from lipid to glycogen-based 
energy metabolism (as well as upregulation of the repro- 
ductive system at the end of the dry season) and whole- 
body fat percentage does not include information on 
these physiological processes. Lastly, the fitness conse- 
quences of (slight) deviation from the phenotypes 
induced by the developmental phenotypic plasticity may 
be higher in the wet season compared to the dry season 
for all or some of the traits. Genetic variation may thus 
be masked and this could involve the activity of molecu- 
lar chaperones that moderate phenotypic expression of 
some genetic variants under wet season conditions more 
than under the dry season conditions (reviewed in Ruth- 
erford 2003; Schlichting 2008). 

The role of phenotypic plasticity was the main reason 
for studying genetic variation for life-history traits in this 
species. It seems evident that in B. anynana, the striking 
plasticity as exhibited in the polyphenism contributes 
much more to the phenotypic variance across seasons 
than genetic variation for the traits (see Fig. 2). In a pre- 
vious study, we have shown that thorax ratio, RMR, and 
to a lesser extent fat content, but not pupal mass and 
developmental time, correlate with the hormonal titers 
that underpin the adaptive phenotypic plasticity (Oostra 
et al. 2011). It is notable that the latter two traits do not 
show GEI in this study. This is also true for RMR, but 
this is likely to be due to the low values of genetic varia- 
tion for this trait in general. We thus find significant GEI 
for the traits intimately related to the adaptive phenotypic 
plasticity, thorax ratio, and fat content. For these traits, 
within one season the correlation between control and 
food-limited condition was much higher in the dry season 
compared to the wet season, which may be related to the 
lower genetic variance estimates in the latter. More inter- 
estingly, cross-environmental correlations between seasons 
within one food treatment are lower for the control ver- 
sus the food limitation condition. This may indicate that 
the processes that contribute to food stress responses are 
similar across the seasons and in a sense override the 
processes involved in the polyphenism. 

As stated before, the generally positive values of the 
cross-environmental genetic correlations suggest that 
expression of similar genes contributes to the genetic vari- 
ance in both seasons and conditions. This, together with 
the observed patterns in the cross-environmental correla- 
tions, indicates that gene expression studies in the differ- 
ent environments and within families will allow for the 
identification of the genes and processes that underpin 
the plasticity and stress responses. This would facilitate 
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the dissection of the composite measure of heritability 
into its functional components. 

Conclusion 

Environmental conditions clearly modulate phenotypic 
variation of life-history traits in B. anynana in addition to 
genetic variances. Interestingly, the extent of these effects 
depends on the trait under consideration, as was demon- 
strated here by very trait-specific patterns in variance com- 
ponents under both seasons and for the two feeding 
conditions. In three of the five traits examined here, the 
genetic variances were little affected by changes in season 
or feeding conditions. For the two traits that show evidence 
of genotype-environment interactions (thorax ratio and fat 
percentage), the main pattern is one in which heritability 
increases under conditions that mark the more adverse 
(dry) season with generally positive genetic correlations 
between seasons. This is in contrast to what is expected on 
the basis of the hypothesis that poor environmental condi- 
tions directly constrain the expression of genetic potential 
and the hypothesis that low cross-environmental genetic 
correlations underlie changes in additive genetic variance 
between seasons (Charmantier and Garant 2005). Taken 
together, this pattern is consistent with low across environ- 
mental variation in important traits including developmen- 
tal time and pupal mass. There is clear seasonal plasticity in 
traits which may reflect differential allocation of resources 
(thorax ratio and fat percentage), possibly reflecting 
season-specific challenges on these allocation decisions and 
an important role for phenotypic plasticity in setting the 
trait values. Experimentally induced food limitation had 
relatively little effect at the genetic level, although some 
traits showed strong phenotypic responses. The generally 
high genetic covariances between the seasons and the food 
limitation manipulations suggest that differences in addi- 
tive genetic variances are mostly due to differential regula- 
tion of the same set of genes; a hypothesis we will be 
testing this in our future research. Our data and analyses 
show that studying the environmental dependency of 
genetic variation and cross-environmental correlations in a 
species with a legacy of selection for adaptive phenotypic 
plasticity has clear advantages as it allows the interpretation 
of the patterns in the light of the species ecology and allows 
for the formulation of testable hypotheses. 
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